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Abstract 

. Dynamic spectrum access in cognitive radio networks can greatly improve the spectrum utilization 

efficiency. Nevertheless, interference may be introduced to the Primary User (PU) when the Secondary 
Users (SUs) dynamically utilize the PU's licensed channels. If the SUs can be synchronous with the PU's 
0^ ', time slots, the interference is mainly due to their imperfect spectrum sensing of the primary channel. 

Q\ . However, if the SUs have no knowledge about the PU's exact communication mechanism, additional 

interference may occur In this paper, we propose a dynamic spectrum access protocol for the SUs 
confronting with unknown primary behavior and study the interference caused by their dynamic access. 
Through analyzing the SUs' dynamic behavior in the primary channel which is modeled as an ON- 
OFF process, we prove that the SUs' communication behavior is a renewal process. Based on the 
■ Renewal Theory, we quantify the interference caused by the SUs and derive the corresponding close-form 

' expressions. With the interference analysis, we study how to optimize the SUs' performance under the 

constraints of the PU's communication quality of service (QoS) and the secondary network's stability. 
Finally, simulation results are shown to verify the effectiveness of our analysis. 

Index Terms 

Cognitive radio, dynamic spectrum access, interference analysis, renewal theory. 
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I. Introduction 

Cognitive radio is considered as an effective approach to mitigate the problem of crowded electro- 
magnetic radio spectrums. Compared with static spectrum allocation, dynamic spectrum access (DSA) 
technology can greatly enhance the utilization efficiency of the existing spectrum resources HI. In DSA, 
devices with cognitive capability can dynamically access the licensed spectrum in an opportunistic way, 
under the condition that the interference to the communication activities in the licensed spectrum is 
minimized [2J. Such cognitive devices are called as Secondary Users (SUs), while the licensed users as 
Primary Users (PUs) and the available spectrum resource for the SUs is referred to as "spectrum hole". 

One of the most important issues in the DSA technology is to control the SUs' adverse interference 
to the normal communication activities of the PU in licensed bands One way is to strictly prevent 
the SUs from interfering the PU in both time domain and frequency domain lH, and the other approach 
is to allow interference but minimizing the interference effect to the PU |5|. For the latter approach, 
the key problem is to model and analyze the interference caused by the SUs to reveal the quantitative 
impacts on the PU. Most of the existing works on interference modeling can be categorized into two 
classes: spatial interference model and accumulated interference model. The spatial interference model 
is to study how the interference caused by the SUs varies with their spatial positions ISHTUll, while the 
accumulated interference model focuses on analyzing the accumulated interference power of the SUs at 
primary receiver through adopting different channel fading models such as what discussed in [91 ifTOl with 
exponential path loss, and in ifTTl |[T2l |[T3l with both exponential path loss and log-normal shadowing. 
Moreover, in |[T4l|[T5ll[T6lll[T7ll . the SUs are modeled as separate queuing systems, where the interference 
and interactions among these queues are analyzed to satisfy the stability region. 

However, most traditional interference analysis approaches are based on aggregating the SUs' 
transmission power with different path fading coefficients, regardless the communication behaviors of 
the PU and the SUs. In this paper, we will study the interference through analyzing the relationship 
between the SUs' dynamic access and the states of the primary channel in the MAC layer. Especially, 
we will focus on the situation when the SUs are confronted with unknown primary behavior. If the SUs 
have the perfect knowledge of the PU's communication mechanism, the interference is mainly from the 
imperfect sensing which has been well studied ifTSl . To simplify the analysis and give more insights 
into the interference analysis, we assumed perfect spectrum sensing in this paper. We show that the SUs' 
dynamic behavior in the primary channel is a renewal process and quantify the corresponding interference 
caused by the SUs' behavior based on the Renewal Theory [19]. There are some works using renewal 
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theory for cognitive radio networks. In f20'|, the primary channel was modeled as an ON-OFF renewal 
process to study how to efficiently discover spectrum holes through dynamically adjusting the SUs' 
sensing period. As the extension works of |[20l , Xue et. al. designed a periodical MAC protocol for 
the SUs in |[2TI . while Tang and Chew analyzed the periodical sensing errors in [22|. In Il23lll24l . the 
authors discussed how to efficiently perform channel access and switch according to the residual time 
of the ON-OFF process in the primary channel. Based on the assumption that the primary channel is 
an ON-OFF renewal process, the delay performance of the SUs were analyzed in Il25lll26l . However, all 
these related works have only modeled the PU's behavior in the primary channel as an ON-OFF process. 
In this paper, we further show and study the renewal characteristic of the SUs' communication behavior 
and analyze the interference to the PU when they dynamically access the primary channel. 

The main contributions of this paper are summarized as follows. 

1) We propose a dynamic spectrum access protocol for the SUs under the scenario that they have no 
knowledge about the PU's exact communication mechanism. By treating the SUs as customers and 
the primary channel as the server, our system can be regarded as a queuing system. 

2) Different from the traditional interference analysis which calculates the SUs' aggregated signal 
power at the primary receiver in the physical layer, we introduce a new way to quantify the 
interference caused by the SUs in the MAC layer. This interference quantity represents the 
proportion of the periods when the PUs' communication are interfered by the SUs' dynamic access. 

3) We prove that the SUs' communication behavior in the primary channel is a renewal process and 
derive the close-form expressions for the interference quantity using the Renewal Theory. 

4) To guarantee the PUs' communication quality of service (QoS) and maintain the stability of 
the secondary network, we formulate the problem of controlling the SUs' dynamic access as an 
optimization problem, where the objective function is to maximize the SUs' average data rate with 
two constraints: the PU's average data rate should not be lower than a pre-determined threshold 
and the SUs' arrival interval and transmission time should satisfy the stability condition. 

The rest of this paper is organized as follows. Firstly, our system model is described in Section JI] 
Then, we present the proposed dynamic spectrum access protocol for the SUs in Section |llll We derive 
the close-form expressions for the interference quantity of two different scenarios in Section |IV] and |V] 
respectively. In Section IVll we discuss how to optimize the SUs' communication performance according 
to the interference analysis. Finally, simulation results are shown in Section IVlIl and conclusion is drawn 
in Section lVlIIl 
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II. System Model 

A. Network Entity 

As shown in Fig.[Tl we consider a cognitive radio network with one PU and M SUs operating on one 
primary channel. The PU has priority to occupy the channel at any time, while the SUs are allowed to 
temporarily access the channel under the condition that the PU's communication QoS is guaranteed. An 
important feature of our system is that the communication mechanism in the primary network is private, 
i.e., the SUs have no knowledge when the PU's communication will arrive. 

For the secondary network, M SUs form a group under the management of one coordinator. The 
coordinator is in charge of observing the PU's behavior, deciding the availability of the primary channel, 
coordinating and controlling the SUs' dynamic access. There is a control channel for command exchange 
between the ordinary SUs and the coordinator. The SUs need to opportunistically access the primary 
channel to acquire more bandwidth for high data rate transmission, e.g., multimedia transmission. 
Considering that the ordinary SUs are usually small-size and power-limit mobile terminals, spectrum 
sensing is only performed by the coordinator. Meanwhile, we assume that all the ordinary SUs are 
half-duplex, which means that they cannot simultaneously transmit and receive data packet. 

B. Primary Channel State Model 

Since the SUs have no idea about the exact communication mechanism of the primary network and 
hence cannot be synchronous with the PU, there is no concept of "time slot" in the primary channel 
from the SUs' points of view. Instead, the primary channel just alternatively switches between ON state 
and OFF state, as shown in Fig.|2] The ON state means the channel is being occupied by the PU, while 
the OFF state is the "spectrum hole" which can be freely occupied by the SUs. 

We model the length of the ON state and OFF state by two random variables Tqn and Tqff 
respectively. According to different types of the primary services (e.g., digital TV broadcasting or cellular 
communication), Tqn and Tqff statistically satisfy different distributions. In this paper, we assume that 
Ton and Tqff are independent and satisfy exponential distributions with parameter Ai and Aq respectively, 
denoted by fof^{t) and /offC^) as follows 

Ton ~ /on (t) = i-e-*/^^ 

(1) 

^ ToFF~/oFF(t) = i;^e-*/^°. 
In such a case, the expected lengths of the ON state and OFF state are Ai and Aq accordingly. These 
two important parameters Ai and Aq can be effectively estimated by a maximum likelihood estimator 
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ll20l . Such an ON-OFF behavior of the PU is a combination of two Poisson process, which is a renewal 
process |[T9]| . The renewal interval is Tp — Ton + Toff SLXid the distribution of Tp, denoted by fp(t), is 

fp{t) = /0N(t) * /oFF(t), (2) 

where the symbol "*" represents the convolution operation. 

III. Secondary Users' Dynamic Spectrum Access Protocol 

In this section, we will design and analyze the SUs' communication behavior including how the SUs 
dynamically access the primary channel and how the coordinator manages the group of SUs. Based on 
the behavior analysis, we can further study the interference caused by the SUs' access. 

A. Dynamic Spectrum Access Protocol 

In our protocol, the SUs who want to transmit data must first inform the coordinator with a request 
command, which can also be listened by the corresponding receiver. The coordinator sequentially responds 
to the requesting SUs by a confirmation command according to the First-In-First-Out (FIFO) rule. The 
SU who has received the confirmation should immediately transmit data with time Tt over the primary 
channel. During the whole process, all the spectrum sensing and channel estimation works are done by 
the coordinator simultaneously. The proposed dynamic access protocol for both ordinary SUs and the 
coordinator is summarized in Algorithm [T] Considering the existence of malicious SUs who may keep 
requesting for packet transmission, we restrict each SU's overall times of requesting within a constant 
period of time. Once the coordinator discovered that one SU's requesting frequency exceeds the pre- 
determined threshold, it will reject this SU's request within a punishing period. 

B. Queuing Model 

According to the proposed access protocol, the secondary network can be modeled as a queueing 
system as shown in Fig.[3] We assume that the requests from all SUs arrive by a Poisson process at 
the coordinator with rate A^^. In such a case, the arrival intervals of SUs' requests at the coordinator, 
denoted by Tg, satisfies the exponential distribution with expectation As, i.e., Tg ~ fs{t) = j-e~^^'^'. 

In this queuing system, the coordinator's buffer only records the sequence of the SUs' request, instead 
of the specific data packets. The packets are stored in each SU's own data memory, which is considered 
as infinite length. In such a case, we can also regard the buffer in the coordinator as infinite length. For 
the service time of each SU, it is the sum of the transmission time Tj and the waiting time if Tt ends 
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Algorithm 1 SUs' Dynamic Spectrum Access Protocol. 

/. For the ordinary SU 

1: if A SU has one packet to transmit tlien 

2: • Send a request command to the coordinator through the secondary control channel 
3: wliUe No confirmation from the coordinator do 
4: • Store the packet in its memory and wait 

5: end wliile 

6: • Transmit its packet after confirmation 

7: end if 

//. For the coordinator 

1: • Estimate the primary channel's parameters Aq and Ai 

2: if A request command is received tlien 

3: • Register the SU's request on the waiting list according to the FIFO rule 

4: end if 

5: wliile The primary channel is in the OFF state do 

6: • Response to the requesting SU on the top of the waiting list with a confirmation command 
7: • Wait for transmission time Tt 

8: end wliile 

9: • Keep sensing the primary channel 



in the ON state as show in Fig.[3] In our model, the time consumed by command exchange between 
ordinary SUs and the coordinator is not taken into account, since it is negligible compared to Tt, Aq and 
Ai- Based on this queuing model, we can analyze the interference caused by the SUs' dynamic access. 

C. Interference Quantity 

If the SUs have the perfect knowledge of communication scheme in the primary network, e.g. the 
primary channel is slotted and all SUs can be synchronous with the PU, then the SUs can immediately 
vacate the occupied channel by the end of the slot. In such a case, the potential interference only comes 
from the SUs' imperfect spectrum sensing. However, when an SU is confronted with unknown primary 
behavior, additional interference will appear since the SU may fail to discover the PU's recurrence 
when it is transmitting packet in the primary channel, as shown by the shaded regions in Fig.[3] The 
essential reason is that the half-duplex SUs cannot receive any command from the coordinator during 
data transmission or receiving. Therefore, the interference under such a scenario is mainly due to the 
SUs' failure of discovering the PU's recurrence during their access time. 

In most of the existing works ||5Tl- |[T7l . interference to the PU was usually measured as the quantity of 
SUs' signal power at primary receiver in the physical layer. In this paper, we will measure the interference 
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quantity based on communication behaviors of the PU and SUs in the MAC layer. The shaded regions 
in Fig.|3] indicate the interference periods in the ON state of the primary channel. In order to illustrate 
the impacts of these interference periods on the PU, we define the interference quantity Qi as follows. 

Definition 1: The interference quantity Qi is the proportion of accumulated interference periods to the 
length of all ON states in the primary channel within a long period time, which can be written by 



In the following sections, we will derive the close form of Qi in two different scenarios listed below. 

• Q/j : SUs with arrival interval = 0. 

• Qi^ : SUs with constant arrival interval As ^ 0. 



In this section, we will discuss the interference to the PU when the average arrival interval of all SUs' 
requests As = 0, i.e., the arrival rate Aj^ = +oo. In the practical scenario, Aj^ = +oo is corresponding 
to the situation when each SU has plenty of packets to transmit, resulting in an extremely high arrival 
rate of all SUs' requests at the coordinator, i.e., A~^ — > +oo. In such a case, the coordinator's buffer is 
non-empty all the time, which means the SUs always want to transmit packets in the primary channel. 
Such a scenario is the worst case for the PU since the maximum interference from the SUs is considered. 

A. SUs' Communication Behavior Analysis 

Since As = means the coordinator always has requests in its buffer, the SUs are either transmitting 
one packet or waiting for the OFF state. As shown in Fig.|4]-(a), the SUs' behavior dynamically switches 
between transmitting one packet and waiting for the OFF state. The waiting time, denoted by T^, will 
appear if the previous transmission ends in the ON state, and the value of is determined by the 
length of the remaining time in the primary channel's ON state. As we discussed in Section IIII-CI the 
interference to the PU only occurs during the SUs' transmission time Tf. Therefore, the interference 
quantity is determined by the occurrence probability of Tt. In the following, we will analyze the SUs' 
communication behavior based on Renewal Theory. 

Theorem 1: When the SUs' transmission requests arrive by Poisson process with average arrival interval 
As = 0, the SUs' communication behavior is a renewal process in the primary channel. 

Proof: As shown in Fig.|4]-(a), we use Ti, to denote the interval of two adjacent transmission 
beginnings, i.e., Tj, = Tt + Tw According to the Renewal Theory |19], the SUs' communication behavior 



Qi = lim 



^ Interference periods 

T 



(3) 




IV. Interference Caused by SUs with Zero Arrival Interval 



November 22, 2011 



DRAFT 



8 



is a renewal process if and only if T^i, Th2, ... is a sequence of positive independent identically distributed 
(i.i.d) random variables. Since the packet transmission time Tt is a fixed constant, Theorem 1 will hold 
as long as we can prove that all Tyji, T^2 ... are i.i.d. 

On one hand, if Tu ends in the OFF state, the following waiting time T^i will be 0, such as T^2 and 
T^4 in Fig.|4]-(a). On the other hand, if Tu ends in the ON state, the length of T^j will depend on when 
this ON state terminates, which can be specifically illustrated in Fig.|4l-(b). In the second case, according 
to the Renewal Theory ||T9ll . is equivalent to the forward recurrence time of the ON state, Tqn, the 
distribution of which is only related to that of the ON state. Thus, we can summarize T^i as follows 

Tti ends in the OFF state, 

ToNi Tti ends in the ON state. 

From (lU, it can be seen that all T^jS are identically distributed. Meanwhile, since each T^i is only 
determined by the corresponding Tu and Toni, all T^jS are independent with each other. Thus, the sequence 
of the waiting time T^i, T^2 ... are i.i.d, which means all Tf,!, T52 . . . are also i.i.d. Therefore, the SUs' 
communication behavior is a renewal process. ■ 



B. Interference Quantity Analysis 

In order to analyze the interference during the SUs' one packet transmission time Tt, we first introduce 
a new function, I{t), defined as follows. 

Definition 2: I{t) is the expected accumulated interference to the PU within a period of time t, where 
t has two special characteristics listed as follows 

• period t always begins at the OFF state of the primary channel, 

• during t, the SUs keep transmitting packets in the primary channel. 

According to Definition 1, Definition 2 and Theorem 1, the interference quantity Qj^ (when = 0) 
can be calculated by 

_ IjTt) 

where /(Tj) is the expected interference generated during the SUs' transmission time Tt, lE(T^) is the 
expectation of SUs' waiting time T^, during which the primary channel is always in the ON state and no 
interference from the SUs occurs. In the following, we will derive the close-form expressions for I{Tt) 
and E(T^) respectively. 
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1} Expected interference I{Tt): According to Definition 2, I{t) is the expected length of all ON states 
within a period of time t, given that t begins at the OFF state. According to the Renewal Theory lfT9l . the 
PU's ON-OFF behavior is a renewal process. Therefore, we can derive I{t) through solving the renewal 
equation Q according to the following Theorem 2. 

Theorem 2: I{t) satisfies the renewal equation as follows 

ct 



lit) = XiFpit) + [ lit- w)fpiw)dw, 
Jo 



(6) 



where fpit) is the p.d.f of the PU's renewal interval given in Q and -Pp(t) is the corresponding cumulative 
distribution function {c.d.f). 

Proof: Let X denote the first OFF state and Y denote the first ON state, as shown in Fig.|5l Thus, 
we can write the recursive expression of function /(i) as follows 

t<X, 

lit) = \ t-X X <t <X + Y, (7) 

Y + lit - X - Y) X + Y<t, 

where X - /opp(x) = ^e'^'/Ao ^nd Y ~ foniv) = 176"^/^^ 

Since X and Y are independent, their joint distribution fxrix^y) = /off (a;) /on (y)- In such a case, 
lit) can be re-written as follows 

^^^^^11 ~ + JJ [y + ^i* - ^ - y)]fxYix,y)dxdy, 

x<t<x+y 2;+y<t 

= J it- x)foj:j:ix)dx+ J J lit- X- y)foPfix)foNiy)dxdy -J J (t-x - y)foffix)fotiiy)dxdy, 
= hit) + hit) - hit), (8) 

where hit), hit) and hit) represent those three terms in the second equality, respectively. By taking 
Laplace transforms on the both sides of ([D, we have 

I(s) = Ii(s) + 12(5) + 13(5), (9) 

where Ii(s), l2is), ^sis) are the Laplace transforms of hit), hit), hit), respectively. 
According to the expression of hit) in (O, we have 

rt 



hit) = it- x)fo^ix)dx = t* fowit). 
Jo 



(10) 
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Thus, the Laplace transform of Ii{t), is 



= ^IFoffCs), (11) 

where Foff(s) = is the Laplace transform of /off(0- 

With the expression of hit) in ([8]l, we have 

hit) = jj lit-x- y)U,ix)fo4y)dxdy = /(t) * ^t) * fo^t) = Ht) * fpit), (12) 

x+y<t 

where the last step is according to Q. Thus, the Laplace transform of hit), his) is 

his) = lis)¥pis), (13) 

where and Fp(s) = ^^^^^^^•''^^^^^j^^ are Laplace transforms of and fpit), respectively. 

Similar to ([T2l ). we can re-written hit) as hit) = t * fpit). Thus, the Laplace transform of hit), 
his) is 

his) = ^Fp(s). (14) 

By substituting ([IT]), and (O into we have 

1 1 F Cs) 

lis) = ^Foff(s) + I(s)Fp(s) - -2Fp(s) = + I(s)Fp(s). (15) 

Then by taking the inverse Laplace transform on the both sides of ( fTSl ). we have 

lit) = Xi [ fpiw)dw + [ lit- w)fpiw)dw = XiFpit) + [ lit- w)fpiw)dw. (16) 
Jo Jo Jo 

This completes the proof of the theorem. ■ 
Theorem 2 illustrates the renewal characteristic of lit). By substituting Fp(s) = into 
(ITSl) . the Laplace transform of can be calculated by 

X,¥pis) Ai 
s(l-Fp(s)) s^(AoAis + Ao + Ai) 



Then, by taking inverse Laplace transform on (fTTl ). we can obtain the close-form expression for l{t) as 

Ao + Ai (Ao + Ai)^v / 



November 22, 2011 



DRAFT 



11 



2) Expected waiting time E(Tu,).- The definition of waiting time T^, has been given in Q in the proof 
of Theorem 1. To compute the expected waiting time, we introduce a new function defined as follows. 

Definition 3: PonC^) is the average probability that a period of time t begins at the OFF state and ends 
at the ON state. 

According to Definition 3 and (|4]|, the SUs' average waiting time ¥,{Tw) can be written as follows 

E(T^) = i^oN(Tt)-E(fo,). (19) 

In the following, we will derive the close-form expressions for i-*oN(7t) and E(Ton), respectively. 

Similar to the analysis of I{t) in Section |IV-B1[ Po^it) can also be obtained through solving the 
renewal equation (l20l ) according to the following Theorem 3. 

Theorem 3: Po^{t) satisfies the renewal equation as follows 

PoN(t) = Ai/p(t) + / Pon{t-w)fp{w)dw. (20) 

^0 

Proof: Similar to I{t) in ([7]), the recursive expression of PoN(i) can be written by 

t<X, 

PoN{t) = \ I X <t < X + Y, (21) 

Pof,{t-X-Y) X + Y<t. 
where X and Y are same with those in In such a case, PoN(i) can be re-written by 



PoN{t) = J J fxY{x,y)dxdy + J J PoN{t - X - y)fxY{x,y)dxdy, 

x<t<.x+y ^+y'^i 

= Fopp(t)-/oFF(t)*i^oN(t) + PoN(i)*/p(i)- (22) 



By taking Laplace transform on the both sides of (I22b . we have 

Pon(s) = AiFon(s) * Fp(s) + Pon(s) * Fp(s). (23) 



Then, by taking the inverse Laplace transform on (1231 ). we have 



PonW = Ai/p(t) + / Po^{t-w)fp{w)dw. (24) 
Jo 

This completes the proof of the theorem. ■ 
Similar to the solution to renewal equation (O in Section IIV-B 1 [ we can obtain the close-form 
expression of Pon(0 through solving (l24l ) as follows 

^oN(t) = ^^^(l-e ^M. (25) 
Ao + '^l ^ ' 
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The Ton is the forward recurrence time of the primary channel's ON state. Since all ON sates follow 
a Poisson process. According to Renewal Theory |fT9ll . we have 

roN~^e-*/^S E(foN) = Ai. (26) 



By combining (l25l) and (l26l) . the SUs' average waiting time E(T^) can be obtained as follows 



IE(7;.) = T^fl - )• (27) 

'^O + '^l ^ ' 

Finally, by substituting (ITSl) and (l27l ) into (O, we can obtain the quantity of interference Qi^ as follows 

(Ao + Ai)Tt-AoAifl-e 



Qh - 7 Aq + Ai t^n ■ (28) 

(Ao + Ai)Tt + A2(l-e~ -om 



V. Interference Caused by SUs with Non-Zero Arrival Interval 

In this section, we will discuss the case when the SUs' requests arrive by a Poisson process with 
average arrival interval As ^ 0. Under such a scenario, the buffer at the coordinator may be empty 
during some periods of time. Similar to the analysis in Section |lVl we will start with analyzing the SUs' 
communication behavior, and then quantify the interference to the PU. 

A. SUs' Communication Behavior Analysis 

Compared with the SUs' behavior when A^ = 0, another state that may occur when A^ ^ is there is 
no SUs' request in the coordinator's buffer. We call this new state as an idle state of the SUs' behavior, 
while the opposite busy state refers to the scenario when the coordinator's buffer is not empty. The length 
of the idle state and busy state are denoted by Tj and respectively. As shown in Fig.|6j the SUs' 
behavior switches between the idle state and busy state, which is similar to the PU's ON-OFF model. In 
the following, we prove that the SUs' such idle-busy switching is also a renewal process. 

Theorem 4: When the SUs' transmission requests arrive by Poisson process with constant rate A^^, 
the SUs' communication behavior is a renewal process in the primary channel. 

Proof: In Fig.|6j we use Tc to denote one cycle of the SUs' idle and busy state, i.e., Tc = T/ + Tb- 
To prove Theorem 4, we need to show that all cycles Td, Tc2, ... are i.i.d. 

For each idle state, its length T/ = Tg is t\\e. forward recurrence time of the SUs' arrival interval Tg 
defined in Section IIII-BI Since the SUs' requests arrive by Poisson process, T/ ~ ^e~*/^= according 
to the Renewal Theory |[T9l . Therefore, the lengths of all idle states are i.i.d. For each busy state, 

N 

Tb = Yli '^bi as shown in Fig.[6l where is the number of SUs' transmitting-waiting times during the 

2=1 
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busy state. Since all Thi are i.i.d as proved in Theorem 1, Tbi, Tb2, ■ ■ ■ will also be i.i.d if we can prove 
that the N of all busy states are i.i.d. It is obvious that the N of all busy states are independent since 
the SUs' requests arrive by a Poisson process. In the following, we will focus on proving its property of 
identical distribution. 

To obtain the general distribution expression of N, we start from analyzing the cases with = 1, 2, 
3 in Fig.|7J where Ei represents the number of requests waiting in the coordinator's buffer at the end of 
the Ith Tt, i.e., the time right after the transmission of the SUs' Ith packet. Thus, as shown in Fig.|7]-(a), 
the probabiUty P{N = 1) is 

P{N = 1) = P{E, = 0\Eo = 1) = Pio, (29) 

where Pij denotes the probability that the last Tt ends with i requests in the coordinator's buffer and 
current Tt ends with j requests. More specifically, Pij represents the probability that there are j — i + 1 
requests arriving at the coordinator during the period T^ + Tt. Since the SUs' requests arrive by a Poisson 
process with arrival interval Ts, Pij can be calculated by 

'j-i+l j-i+2 

T^ \ = I 

(j-i + l)! 



= P\Y. Ts.<iT^ + Tt)< E = / )^_'',., e-'/>^^P{T^ + Tt = t)dt, 



k=l k=l 



Tt 



oo ((t + Tt)/X, 



(j-i + l)! 



e-^'+''>'^'P{T^ = t)dt, (30) 



where T^^^ is the SUs' A;th arrival interval satisfying the exponential distribution with parameter A^, the 

j-i+l i-i+2 

first equality is because ^ Tg^ and ^ Ts^ satisfy Erlang distribution. 

k=l k=l 

According to (|4]| and (l26l) . the probability distribution of T^, P{T^ = t), can be written as follows 



P(T^ = t) = { 



PoFF(Tt) t = 0, 

(31) 



where Po,,{Tt) = 1 - Po^{Tt) = ^ (l - e Xm'^* 



By substituting (|3T|) into (I30b , we can re-write Pij as 

{j-i + iy. jo+ + Ai 

When = 2, as shown in Fig.|7]-(b), P{N = 2) can be written by 

P{N = 2) = P{E2 = ^,Ei = 1\Eq = 1). (33) 
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According to the queuing theory |[27l . the sequence Ei, E2, E]\f is an embedded Markov process. 
Thus, P{N = 2) can be re-written by 

P{N = 2) = P{Ei = l\Eo = l)P{E2 = 0|^i = 1) = PiiPio. (34) 

When = 3, as shown in Fig.|7]-(c), there are two cases (£'0 = 1, Ei = 2, E2 = 1, E2 = 0) and 
(Eo = 1, Ei = 1, E2 = 1, E2 = 0). Thus, 

P{N = 3) = P12P21P10 + PuPnPio- (35) 

When N = n, El {0 < I < n) should satisfy the following condition 

Eo = l, 1 < Si < n - 1, . . . , l<Ei<n-l, E^-i = 1, En = 0. (36) 

Thus, there are (n — 1)! possible combinations of {Eq, Ei, .. ., Ei, . .. , En-i, En)- We denote each case 
as C(a), where 1 < a < (n — 1)!. For each case, the probability is the product of n terms Pij(C{a), 6), 
where 1 < 6 < n. Thus, P{N = n) can be expressed as follows 

(n-l)! n 

P{N = n)= W Pij {C{a),b). (37) 

a=l b=l 

From (I37] ). we can see that N of all busy states are identical distributed, and hence i.i.d. 

Up to now, we have come to the conclusion that Tj of all idle states are i.i.d, as well as Tb of all 
busy states. Since Tj and Tb are independent with each other, the sequence of all cycles' lengths Td, 
Tc2, ... are i.i.d. Therefore, the SUs' communication behavior is a renewal process. ■ 

B. Interference Quantity Analysis 

According to Definition 1 and Theorem 4, the interference quantity Qi^ can be calculated by 

Qh=I^B-Qh, (38) 

where hb = e{Ti)+e{Tb) '■^^ occurrence probability of the SUs' busy state. 

Our system can be treated as an M/G/1 queuing system, where the customers are the SUs' data packets 
and the server is the primary channel. The service time S of one SU is the sum of its transmission time Tt 
and the waiting time of the next SU T^. In such a case, the expected service time is E(S') = + E(T^). 
According to the queuing theory |[27l . the load of the server is p = ¥,{S)/X, where A is the average 
arrival interval of the customers. By Little's law |27], p is equivalent to the expected number of customers 
in the server. In our system, there can be at most one customer (SUs' one packet) in the server, which 
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means the expected number of customers is equal to the probability that there is a customer in the server. 
Therefore, p is equal to the proportion of time that the coordinator is busy, i.e.. 

Thus, combining (|28]) . (f38l) and (f39l) . the close-form expression of Qj^ can be obtained as follows 



(Ao + \i)Tt - AoAi (l - e Xm'^') 



= IJmTm) 

VI. Optimizing Secondary Users' Communication Performance 

In this section, we will discuss how to optimize the SUs' communication performance while maintaining 
the PU's communication QoS and the stability of the secondary network. In our system, the SUs' 
communication performance is directly dependent on the expected arrival interval of their packets aJi] and 
the length of the transmission time Tt. These two important parameters should be appropriately chosen 
so as to minimize the interference caused by the SUs' dynamic access and also to maintain a stable 
secondary network. 

We consider two constraints for optimizing the SUs' A^ and Tt as follows 

• the PU's average data rate should be at least Rp, 

• the stability condition of the secondary network should be satisfied. 

In the following, we will first derive the expressions for these two constraints based on the analysis in 
Section |IV] and |Vl Then we formulate the problem of finding the optimal A* and as an optimization 
problem to maximize the SUs' average data rate. 

A. The Constraints 

1 ) PU's Average Data Rate: If there is no interference from the SUs, the PU's instantaneous rate is 

log(l + SNRp), where SNRp denotes the Signal-to-Noise Ratio of primary signal at the PU's receiver. On 

/ SNR \ 

the other hand, if the interference occurs, the PU's instantaneous rate will be \ogyi + jj^j^ ^^^''^ 
INRp is the Interference-to-Noise Ratio of secondary signal received by the PU. According to Definition 
1, Qj^ represents the ratio of the interference periods to the PU's overall communication time. Thus, the 
PU's average data rate Rp can be calculated by 

Rp={l- Q/J • log(l + SNRp) + Qj, ■ log(^l + • (41) 

'To evaluate the stability condition, we only consider the scenario when Aa 7^ 0. 
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2) SUs' Stability Condition: In our system, the secondary network and the primary channel can be 
modeled as a single-server queuing system. According to the queuing theory |[27l . the stability condition 
for a single-server queue with Poisson arrivals is that the load of the server should satisfy p < I 1281 . 
In our system, we have 

p = / ^ < 1. (42) 
In such a case, the SUs' stability condition function, S{Tt,\s), can be written as follows 

A? / hn+±XT \ 

S{Tu\s) = \s-Tt-—^(l-e- >0. (43) 

Aq + Ai V / 

B. Objective Function: SUs' Average Data Rate 

If a SU encounters the PU's recurrence, i.e., the ON state of the primary channel, during its transmission 
time Tt, its communication is also interfered by the PU's signal. In such a case, the SU's instantaneous 
rate is log^l + j^^^p^) ' ^^^^^ SNRg is the SU's Signal-to-Noise Ratio and INR^ is the Interference-to- 
Noise Ratio of primary signal received by the SU. According to Theorem 1 and Theorem 4, the occurrence 
probability of such a phenomenon is P'B j'^^^^x ) ~ ^T^- other hand, if no PU appears during 

the SU's transmission, its instantaneous rate will be log(l + SNR^) and the corresponding occurrence 
probability is PB t'+e{tJ) = aT ■ ^hus, the SU's average data rate Rs is 

Tt - liTt) ( \ I(Tt) ( SNR^ \ 

C. Optimizing SUs' Communication Performance 

Based on the analysis of constraints and objective function, the problem of finding optimal and A* 
for the SUs can be formulated as follows 

, Tt-IiTt) ( \ I(Tt) ( SNR^ 

max R,{Tu A. = * / • log 1 + SNR^ + ^ • log 1 + — — ^ ^ 

{TuK) ^ ' \s V ) \s \ INRs + 1^ 

s.t. R^{Tu Xs) = (1 - Q/.) • log(l + SNRp) + Qi^ ■ log(^l + 5;^^^) > R^, (45) 



At / Ao + Al 7^ \ 

S{Tt,Xs) = Xs-Tt-j-^^{l-e — ^'j>0. 



Theorem 5: The SUs' average data rate Rs{Tt, Xs) is a strictly increasing function in terms of the their 
transmission time Tt and a strictly decreasing function in terms of their average arrival interval A^, i.e., 

£<«■ 
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The PU's average data rate Rp{Tt,Xs) is a strictly decreasing function in terms of Tt and a strictly 
increasing function in terms of A^, i.e., 

dRn dRr> 

The stability condition function S{Tt, As) is a strictly decreasing function in terms of Tj and a strictly 
increasing function in terms of A^, i.e., 

dS dS 

m^"- w.^"- '''' 

Proof: For simplification, we use Rso to express log^l+SNRg^ and Rgi to express logfl+f^^^ 
According to (l44l ) and (fTSl) . and can be calculated as follows 

dRs Rso Rso — Rsi dI{Tt 



INR,+ir 



dT, X, X, dTt 



As(Ao + Ai) 
dR, 1 



^ ^ XoRsQ + XiRsi + Xi{Rso-Rsi)(l-e (49) 



5A, A2 



Tt-I{Tt)]RsQ + I{Tt)Rsi]. (50) 



Since fi^o > -Rsi, e "oa/^' < 1, and Tt > I{Tt), we have 



dRs dRs 

Similarly, we use RpQ to express log^l + SNRp^ and R^i to express log^l + j^^^^^ ■ According 



to (|4TI ). and can be calculated as follows 



(i?sO — -Rsl), -^T— = — (-^^0 ~ -^gl)- (52) 



dRp _ dQi^ dRp dQi^ 



dTt dTt ' '''' dXs dX 

According to ( |40l ). we have 



>0. ^<0. (53) 



dTt Xs ' dX 

Thus, combining (l52l ) and (1531 ). we have 



According to ( |43l) . ^ and can be calculated as follows 

This completes the proof of the theorem. ■ 
From Theorem 5, we can see that the objective function and the constraints are all monotonous functions 
in terms of Tt and A^. Thus, the solution to the optimization problem (|46l ) can be found using gradient 
descent method ||29l. 
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VII. Simulation Results 

In this section, we conduct simulations to verify the effectiveness of our analysis. The parameters of 
primary ON-OFF channel are set to be Aq = 2.6s and Ai = 3.6s. According to Fig.[3j we build a queuing 
system using Matlab to simulate the PU's and SUs' behaviors. 

A. Interference Quantity Qj 

In Fig.[8] and Fig.|9j we illustrate the theoretic and simulated results of Qj^ and Qj^, respectively. 
The theoretic Qj^ and Qj,^ are computed according to (1281 ) and (l40l ) with different values of the SUs' 
transmission time Tf. The average arrival interval of the SUs' packets A^ is set to be 1.3s when calculating 
theoretic Qj^- For the simulated results, once the interference occurs, we calculate and record the ratio 
of the accumulated interference periods to the accumulated periods of the ON states. 

From Fig.[8]and Fig.|9j we can see that all the simulated results of Qj^ and Qj^ eventually converge to 
the corresponding theoretic results after some fluctuations at the beginning, which means that the close- 
form expressions in (|28] ) and ( |40l ) are correct and can be used to calculate the interference caused by the 
SUs in the practical cognitive radio system. Moreover, we can also see that the interference increases as 
the SUs' transmission time Tj increases. Such a phenomenon is because the interference to the PU can 
only occur during Tt and the increase of Tt enlarges the occurrence probability of Tt. Finally, we find 
that due to the existence of the idle state when As / 0, Qj^ is less than Qj^ under the same condition. 

B. Stability of The Secondary Network 

Since we have modeled the secondary network as a queuing system shown in Fig.|3j the stability 
of the network is reflected by the status of the coordinator's buffer. A stable network means that the 
requests waiting in the coordinator's buffer do not explode as time goes to infinite, while the requests in 
the buffer of an unstable network will eventually go to infinite. In Section IVI-A2I we have shown the 
stability condition of the secondary network in (1431 ). On one hand, if the SUs' access time Tt is given 
in advance, the SUs' minimal average arrival interval As can be computed by (1431 ). On the other hand, 
if As is given, the maximal Tt can be obtained to restrict the SUs' transmission time. 

In this simulation, we set Tt = 0.6s, and thus As should be larger than 1.25s to ensure the SUs' 
stability according to (|43] ). In Fig.|9j we show the queuing length, i.e., the number of requests in the 
coordinator's buffer, versus the time. The black lines shows the queuing length of a stable network, in 
which As = 1.3s is larger than the threshold 1.25s. It can be seen that the requests dynamically vary 
between and 60. However, if we set As = 1.2s, which is smaller than the lower limit, from Fig.|9j we 
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can see that the queuing length will finally go to finite, which represents an unstable network. Therefore, 
the stability condition in (1431) should be satisfied to maintain a stable secondary network. 

C. PU's and SUs' Average Data Rate 

The simulation results of the PU's average data rate Rp versus the SUs' transmission time Tt and 
arrival interval As are shown in Fig.[TT] where we set SNRp = SNR^ = 5db and INRp = INR^ = 3db. We 
can see that Rp is a decreasing function in terms of Tt given a certain Xg, and an increasing function in 
terms of for any fixed Tt, which is in accordance with Theorem 5. Such a phenomenon is because an 
increase of Tt or a decrease of A^ will cause more interference to the PU and thus degrade its average 
data rate. In Fig.[T2l we illustrate the simulation results of the SUs' average data rate Rs versus Tt and 
As. Different from Rp, Rs is an increasing function in terms of Tt given a certain As, and a decreasing 
function in terms of As for any fixed Tt, which also verifies the correctness of Theorem 5. 

Suppose that the PU's data rate should be at least 2.0bps/Hz, i.e., Rp = 2.0bps/Hz. Then, according 
to the constraints in (l46l ). Tt should be no larger than the location of those three colored vertical lines 
in Fig.fTT] corresponding to As = 1.3s, 1.5s, 2.0s respectively. For example, when As = 1.3s, the optimal 
Tt should be around 400ms to satisfy both the Rp and stability condition constraints. In such a case, the 
SUs' average data rate can achieve around 0.6bps/Hz according to Fig. [121 For any fixed Rp, the optimal 
values of Tt and A* are determined by the channel parameters Aq and Ai. Therefore, the SUs should 
dynamically adjust their communication behaviors according to the estimated channel parameters. 

VIII. Conclusion 

In this paper, we analyzed the interference caused by the SUs confronted with unknown primary 
behavior. Based on the Renewal Theory, we showed that the SUs' communication behaviors in the 
ON-OFF primary channel is a renewal process and derived the close-form for the interference quantity. 
We further discussed how to optimize the SUs' arrival rate and transmission time to control the level of 
interference to the PU and maintain the stability of the secondary network. Simulation results are shown to 
validate our close-form expressions for the interference quantity. In the practical cognitive radio networks, 
these expressions can be used to evaluate the interference from the SUs when configuring the secondary 
network. In the future work, we will study how to concretely coordinate the primary spectrum sharing 
among multiple SUs. 
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Fig. 1. Network entity. 
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Fig. 2. Illustration of the ON-OFF primary channel state. 




Fig. 3. Illustration of the queuing system, the SUs' dynamic spectrum access and interference to the PU. 
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Fig. 4. Illustration of the SUs' behavior in the primary channel when As = 0. 
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Fig. 7. Illustration of buffer status Ei when = 1, 2, 3. 
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Fig. 9. Interference quantity Qir^ with As = 1.3s. 
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Fig. 10. Queuing length under stable and unstable conditions. 
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